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H ! ABSTRACT 

Context. Due to the star activity, the masses of CoRoT-7b and CoRoT 7c are uncertain. Investigators of the CoRoT team have proposed 
several solutions, all but one of them larger than the initial determinations of 4.8 ± 0.8 MEarth for CoRoT-7b and 8.4 ± 0.9 Meoi± for 
CoRoT 7c. 

Aims. This investigation uses the excellent HARPS radial velocity measurements of CoRoT-7 to re-determine the planet masses and 
to explore techniques able to determine mass and elements of planets discovered around active stars when the relative variation of the 
radial velocity due to the star activity cannot be considered as just noise and can exceed the variation due to the planets. 
Methods. The main technique used here is a self-consistent version of the high-pass filter used by Queloz et al. (2009) in the first mass 
determination of CoRoT-7b and CoRoT-7c. The results are compared to those given by two alternative techniques: (1) The approach 
proposed by Hatzes et al. (2010) using only those nights in which 2 or 3 observations were done; (2) A pure Fourier analysis. In all 
cases, the eccentricities are taken equal to zero as indicated by the study of the tidal evolution of the system; the periods are also kept 
fixed at the values given by Queloz et al. Only the observations done in the time interval BJD 2,454,847 - 873 are used because they 
include many nights with multiple observations; otherwise it is not possible to separate the effects of the rotation fourth harmonic 
(5.91 d=P TOt /4) from the alias of the orbital period of CoRoT-7b (0.853585 d). 
Results. The results of the various approaches are combined to give for the planet masses the values 8.0 + 1.2M Earlll for CoRoT-7b 
and 13.6 ±1.4 ME art h for CoRoT 7c. An estimation of the variation of the radial velocity of the star due to its activity is also given. 
Conclusions. The results obtained with 3 different approaches agree to give masses larger than those in previous determinations. From 
the existing internal structure models they indicate that CoRoT-7b is a much denser super-Earth. The bulk density is 1 1 ± 3.5 g.cirT 3 . 
CoRoT-7b may be rocky with a large iron core. 



CO 

> . 

Key words. star:individual:CoRoT 7 - planetary systems - stanactivity - methods: statistical - techniques: radial velocities - mass 
determination - CoRoT 7b - CoRoT-7c - exoplanets: hot super-Earths 



i~ I ' 1. Introduction variation in the measured radial velocity comes mainly from the 
,— ' , activity of the star whose spots determine the value of the radial 
O . CoRoT-7b was the first super-Earth for which mass and radius ve locity integrated over its disk. This activity follows roughly 
r-H . have been determined. CoRoT-7b and the recently discovered the rotation per i d of the star (23.64 days), but in a very irregu- 
> ! GJ 1214b (Charbonneau et al. (|2009j)) and Kepler-lOb (Batalha lar way . me 150 . day photometric observations done by CoRoT 
• rH . et al. <|2011|» are paradigms for the study of the physics of what show the variation displaying full span in some times, but almost 
[X; ■ exo-Earths, super-Earths and/or mim-Neptunes can be. They set disa pp earing in other s (see Leger et al. (l2009h ). Because of this 
H ; the only real constraints available to models of the formation and activity? CoRoT-7b is also the paradigm of the kind of problems 
y* , evolution of hot telluric planets. For this reason, it is very impor- that may be often found when planets whh Earth _i ike masses ^ 
tant to have good radius and mass determinations. The radius of discover ed. One of the aims of the investigation reported in this 
CoRoT-7b, determined from the transits observed by CoRolU is paper is to use Co RoT-7, for which an excellent set with 109 
10,100 ± 600 km (Bruntt et al. (|2010|)), a value that may be im- HARPS radial velocity measurements exists, to explore tech- 
proved, but whose magnitude is nevertheless definitively estab- niques that may be used when dealing with i ow _ mass planets, 
hshed. The mass, determined from the radial velocity measure- We may wish that future discove ries are done around more quiet 
ments (4.8 + 0.8 Earth masses cf Queloz et al. <|2009D; hereafter starSi but Earth . like planets are too important and we cannot dis- 
QBM), however, has not the same accuracy. This is due to the card any of them because of the central star activity. We hope that 
fact that the 0.91 + O.O3M (see Lanza et al. ([20101)), 1-2 Gyr-old CoRol ; KEPLER and ground based instruments will discover 
G9V star CoRoT-7 (=TYCHO 4799-1733-1) is too active. The new EarthSi super -Earths and mini-Neptunes which, as Kepler- 
10b, may be in orbit around quiet stars. However, as CoRoT-7b, 



1 The CoRoT space mission, launched on December 27th 2006, has they may be found around non-quiet stars and the improvement 
been developed and is operated by CNES, with the partnership of G f the techniques of mass and orbit determination used to study 
Austria, Belgium, Brazil, ESA, Germany and Spain. 
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Fig. 1. Relative radial velocity measurement (dots) and the con- 
tribution of the activity as estimated with a filtering using the first 
three harmonics of the rotational period (solid line). (Adapted 
fromFig.8ofQBM) 

such cases is important. This need is motivating a great deal of 
investigations and the number of p apers dealing with th e mas s 
of CoRoT-7b (QBM; Boisse et al. d2011l) : Hatzes et al. d2010t) ; 
Pont et al. d201 lb ) is increasing. In addition to this interest, we 
have to consider that HARPS, currently the only instrument able 
to make measurements of CoRoT-7 with the required precision, 
is being on demand by a great deal of other targets and we can- 
not foresee when a new series of measurements with the same 
quality of the existing one will become possible. 
The nature and magnitude of the problem of the mass and orbit 
determination of the CoRoT-7 planets can be assessed from fig.Q] 
(adapted from fig. 8 of QBM). It shows the 5 sets of measure- 
ments of the radial velocity of CoRoT-7 obtained with HARPS 
(dots) and the estimated part of the radial velocity due to the star 
activity. The activity shown in fig. Q] (solid line) was obtained 
in QBM by means of a filtering designed to eliminate periodic 
disturbances corresponding to the first three harmonic^ of the 
rotational period. These sets present patterns very different one 
from another showing that the construction of one single model 
for the activity in the whole interval is not possible. It also dis- 
courages the extensive use of Fourier tools (periodograms) over 
the whole set. We have to search for techniques able to separate 
the low-frequency rotation signals (periods 23.64 days and its 
main harmonics), from the higher frequency signal coming from 
the two planets (periods 0.853585 and 3.698 days cf. QBM). The 
difference between the periods of the two components (rotation 
and planets) indicates that, in this case, we may filter the data 
from its low-frequency parts and consider separately the high- 
frequency information. In the case of one continuous signal, or 
at least of a long discrete evenly spaced time series, the problem 
that we have to solve is classical and well known: We should 
construct and use a high-pass filter. However, the simple recipes 
to construct a high-pass filter in the frequency domain cannot 
be used for a series with a limited amount of data and, worse, 
unevenly distributed. Because of the high correlation between 
the absolute values of the Fourier transform (or spectral power) 
at different frequencies, these filters must necessarily be con- 
structed in the time domain. 

One high-pass filter in the time domain was used in QBM in 
the first mass determination of CoRoT-7b and-7c. However, as 
discussed in Section 2 of this paper, the filter then used affected 
the high-frequencies, the importance of which was downsized 
without apparent reason, and set the cut-off below the frequency 

2 We follow the common usage in Physics where N' 1 ' harmonic means 
an oscillation whose frequency is N times the fundamental frequency. 
Thus the fundamental frequency 2n/P Iot is the first harmonic, the sec- 
ond harmonic is the component whose period is P mt /2, and so on. 



of the main alias of the period of CoRoT- 7b, thus affecting also 
the amplitudes corresponding to this planet (cutting out the alias, 
the filter also considerably downsized the signal corresponding 
to the actual frequency of the planet). 

This paper starts with an analysis of the mass determination pub- 
lished in QBM's CoRoT-7c discovery paper (section 2) and then 
proposes a self-consistent agorithm founded on the filtering tech- 
nique used there (section 3), which is applied to the selected set 
of 52 observations made in 27 consecutive nights, between BJD 
2,454,847 and BJD 2,454,873 (section 4) to obtain the masses 
of the planets. The restriction of the analysis to this set of dates 
is due to the fact that, in 10 nights of this set, 3 observations 
were done, spanning about 4 hours between the first and last 
observation in the night. This is not enough to completely de- 
stroy the aliasing due to the almost uniform spacing between 
observations done in consecutive nights, but it allows us to dis- 
tinguish between two solutions with forced p eriods equal to the 
transit period and its alias (see Hatzes et al. d2010h figure 7). It 
is worth stressing the fact that in a series made of observations 
taken always near the same hour in the night, no mathemati- 
cal tool exists able to distinguish between one frequency and its 
aliases. In section 4, the dependence of the results on the filter 
parameters is also discussed. In sections 5 to 7, we present the 
resulting estimate of the star activity and analyze the residuals 
obtained by subtracting the activity from the observed radial ve- 
locities. Alternative techniques are discussed in sections 8 and 
9. textbf The approach discussed in section 8 follows a sug- 
estion by Hatzes et al. (2010) and uses only the observations 
from nights where multiple observations were done. This is of 
particular importance because it allows an analysis independent 
of any explicit hypotheses on the behavior of the star activity. 
These observations are analysed here with the help of a biased 
Monte-Carlo technique allowing confidence intervals to be ob- 
tained. The approach discussed in section 9 is a classical multi- 
period Fourier analysis. It differs from other approaches using a 
Fourier decomposition by the fact that, here, no a priori periods 
are used. The periods of the solution found are those allowing us 
to get the best fit of the observations to a multi-periodic function. 
In section 10, we present some simulation results taking into ac- 
count tidal interactions showing the circularization of the orbits 
and thus justifying the adoption of zero eccentricities for both 
planets. At last, we proceed with the discussion of the results 
and the conclusions. 



2. Analysis of the first mass determination 

The analysis of the mass determination published in the CoRoT- 
7c discovery paper (QBM) is the first step in this study and the 
main point to be considered concerns the filtering properties of 
the procedure used there. Is it equivalent to a high-pass filter? In 
order to know that, we compute the Fourier transforms Q of the 
given data and of the residuals obtained in QBM after subtract- 
ing the activity, respectively, and compare them one to another. 
The result presented in fig. [2] shows that, indeed, the used pro- 
cedure completely filtered the low-frequencies (the transform of 
the filtered series is close to zero for all frequencies below 0.22 
d -1 . However, it also affected the high-frequency components. 
The strong alias of the period of CoRoT-7b, at 5.925d (indicated 



3 Spectra obtained using date-compensated discrete Fourier trans- 
forms (DCDFT; cf Ferraz-Mello j!981h ). DCDFT differs from usual 
Lomb periodograms because they also consider the constant compo- 
nent, whose neglect may affect the height of the p eaks ( see discussion 
on floating-mean periodograms in Cumming et al. dl999h ) 
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Fig. 2. Fourier transform (DCDFT) of the measured radial veloc- 
ities (black) and of the values obtained in QBM after subtracting 
the estimated activity (red). Both transforms are presented us- 
ing the same units. C and B indicate the orbital frequency of 
CoRoT-7c and the 5.925-day alias of the CoRoT-7b frequency, 
respectively 

with B' in fig.|2]i, disappeared. The signal of CoRoT-7c was less 
affected by the filter but it was downsized without apparent rea- 
son, since no low-period terms were included in the filter. The 
main reason for this result is the fact that the filtering was done 
on the raw data, without taking into account the part of the signal 
due to the planets. We may guess that the fact that these periods 
are not at all commensurable with the Coherence Time (window) 
used in the filter (20 days) may have played a role in the deep 
sculpting done by the filter at the planet frequencies and their 
aliases (the average of a periodic function over a time interval 
not commensurable with its period is different of zero.). 
A second factor may have been the fact that fittings are exact 
when the number of dates is smaller (as well known). Thus, near 
the borders of the interval, the estimated activity will be closer 
to the given data than in the middle of the period. The end effect 
seems to be responsible for the fish-like appearance of the fil- 
tered velocities in the largest set (narrow in the extremities and 
wide in the middle; See the grey line and dots in fig. |SJ. One 
could think that such appearance might result from a particular 
beat of the orbital frequency of CoRoT 7c and the nearby alias 
of the orbital frequency of CoRoT-7b, but an a posteriori plot of 
synthetic velocities shows that this is not so. Because of the ac- 
tual phases and periods, the destructive interference of the two 
sine curves (at the actual discrete observation times) does not 
happen in the borders of the interval, but near the middle of it. 
It is worth adding that one unconstrained 3-sinewave analy- 
sis of the filtered radial velocities in QBM, in the period BJD 
2,454,847 - 873, using a genetic algorithm completed with a 
downhill simplex gave as more important periods present in the 
filtered data, 3.495d and 3.963d. The difference in frequency of 
these two periods is 0.0338d~\ which is the inverse of 29.6 d 
(very close to the timespan of the observations used in the anal- 
ysis), clearly showing the interference of the timespan of the ob- 
servations in the considered subset. 



3. A self-consistent high-pass algorithm 

The high-pass filter used in QBM may be shortly described as 
follows. First, we define one time window (the Coherence Time), 
fixed as being a guess on the number of days in which it may 
be reasonable to fit the given harmonic function to the activity. 
Then, we construct N time windows of the chosen size, each 
centered on one of the N dates of the observations and including 



all observations inside the window. In every window, the data 
are Fourier analyzed and represented by the first terms of a har- 
monic series whose fundamental period is the star rotation pe- 
riod (23.64 d) plus a constant. The activity at a given date is es- 
timated as the (weighted) average of the values given to it by the 
harmonic representations of the signal in all subsets including 
the given date. For the sake of clearness, let us add the following 
information: (i) The actual window size is fixed in such a way 
that the ends of the window do not separate observations done in 
the same night (if necessary, the actual window is taken slightly 
larger than the nominal Coherence Time); (ii) Near the borders, 
the subsets are incompletely filled as the windows extend to be- 
yond the considered interval covering some nights where no ob- 
servations were done. 

In this paper, we propose an improvement of this procedure. The 
main change is that, now, a predicted signal corresponding to the 
two planets is subtracted from the radial velocities before the fit- 
ting of the harmonic function. Next, the activity is estimated as 
described above and subtracted from the observed radial veloc- 
ities; the resulting residuals are used to determine the masses. 
The new masses may then be used as first guesses in a new run 
of the algorithm and leads to an improved prediction of the sig- 
nal corresponding to the two planets, and so on. The procedure 
is iterated as many times as necessary up to reach a satisfactory 
convergence. Formally, we may say that the whole procedure de- 
fines a map x n+ \ = F(x n ), (here, x represents the masses) which 
is iterated up to get x n+ i - x n . One problem appearing in the 
actual application of this scheme is a possible slow convergence 
of this map. For this reason, to obtain the results discussed in 
this paper, we have rather used an alternative accelerated map: 
x n+x = F(x„ + A(x n - x„_i)) with A = 1. 

In a self-consistent solution, the values of the masses used, to 
subtract the planets, is equal to that obtained from the filtered se- 
ries. Figure[3]shows the evolution of the masses obtained in five 
chains of determinations using the same high-pass filter (with 
4 harmonics and the same Coherence Time). The longest one 
was obtained iterating the process from a first approximation in 
which the actual observed radial velocities are used without tak- 
ing into account the planets (i.e. starting with the two masses 
equal to zero). The others started from arbitary sets of values 
corresponding to higher masses, represented by diamonds in fig. 
[3] One may note the very slow convergence of the iterative pro- 
cedure. One may also note that the algorithm converges quickly 
to a line that corresponds to a linear relation between the two 
masses in which one of them (mj c ) is almost constant. To have 
different convergence ratios along two orthogonal directions is 
a common feature in maps; in this case it seems to be due to 
the slow separation between the planet CoRoT 7b and the fourth 
rotation harmonic. However, notwithstanding the slow conver- 
gence, the five runs converged to the same point P, represented 
by a star in Fig.[3]The reliability of the map was checked using 
some synthetic sets of data constructed using the first, second 
and fourth rotation harmonics with the amplitudes indicated in 
some Fourier analyses (see section 9), two sine curves of half- 
amplitude 6 m.s~', corresponding to two planets, and a Gaussian 
noise. The results are as follows: (1) When the noise is not in- 
cluded, the map reproduces the two planets exactly as given; 
(2) When noise in the range 1.7-2.5 m.s -1 is added, the re- 
sults for CoRoT-7b fall around the given value within 0.3 m.s -1 . 
However, the results for CoRoT-7c fall systematically ~1 m.s -1 
below the given value. We have taken these results into account 
in estimating the error bars of the final results. 
The analysis of the spectrographic parameters related to the 
star activity (FWHM, Rhk, bisector span) shows the contribu- 
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Fig. 3. Evolution of the masses in five chains of results iteratively 
obtained from five different sets of initial conditions (diamonds) 
using the 4-harmonic high-pass filter and observations done in 
the period BJD 2,454,847 - 873. All of them converge to the 
same point P (star). One of the chains started at (0,0). The figure 
also shows the result of the first iteration in QBM (point Ql) and 
the mass values adopted in that paper. 



tion of the rotation higher harmonics to the observed activity. 
However, they were of little help to decide on the kind of filter- 
ing to be used. The power spectra of these parameters in the pe- 
riod BJD 2,454,847 - 873 are shown in the fig.g] Unfortunately, 
the features of the spectra of the bisector span and FWHM show 
only small bumps at the position of the higher-order harmon- 
ics, but similar bumps were found when using scrambled data 
thus showing that the observed one cannot be distinguished from 
bumps generated by white noise and are thus meaningless. The 
corresponding data seems to be affected by the low brightness 
of the star. The only power spectrum showing significant peaks 
above the minimum level of significance is the power spec- 
trum of the index log where peaks corresponding to higher- 
order harmonics are clearly seen even if some offset due to the 
short time span of the observations used can be noted. A simi- 
lar a nalysis using the photometric observations done by Queloz 
et al. (l2009l) (QBM) showed almost no influence of the 4 th har- 
monic and led them to neglect it in the construction of the high- 
pass filter. We will use both filters with and without the fourth 
harmonic and consider them in the composition leading to the 
conclusions of this paper on the mass of the planets. In addition 
we mention that the power spectra clearly show some higher- 
order harmonics which may affect our results. The consideration 
of them would require new improvements, different of those pre- 
sented in this paper. 

We shall mention that the use of the rotation's fourth harmonic 
raises some critical questions. Indeed its period is one fourth of 
23.64d, that is 5.917 d, and one of the main aliases of the pe- 
riod of 7b is 5.925 d. However, aliases are defined for uniformly 
spaced time series with observations separated by multiples of 
one constant value (e.g. one sidereal day) and cannot be avoided 
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Fig. 4. Power spectra of the spectrographic parameters associ- 
ated with the star activity: bisector span, FWHM and log R^, in 
the period BJD 2,454,847 - 873. The exact location of the main 
rotation harmonics is indicated. The gray rectangle is below the 
minimum significance level. Peaks inside this rectangle can be 
produced by random data. 

as far as the separation between observations is kept unaltered. 
This is a classically known fact and it has been taken into ac- 
count in the scheduling of the spectroscopic observations as soon 
as the observations showed this coincidence. The observations 
done in the period JD 2,454,847 -873 include 52 data obtained 
in 27 consecutive nights and in 10 of these nights 3 observa- 
tions were done covering about 4 hours. This was not sufficient 
to completely eliminate aliasing problems, but power spectra ex- 
tended beyond the nominal Nyquist frequency showed that the 
perfect mirroring of the power spectrum has been avoided ( see 
fig. 7; see fig. 10 of QBM). Several other tests were done. Least 
squares determinations involving simultaneously the photomet- 
ric period 0.853585 d and the fourth rotation harmonic lead to 
correlation values in the range 0.70 - 0.88 depending on the de- 
sign of the experiment done and the observations used. Some of 
these results are worrisome, but the final test is provided by the 
coincidence of the convergence points in chains starting at very 
different mass values (as the test ones shown in fig- [3j. One char- 
acteristic of least-squares procedures involving highly correlated 
parameters is that the results become erratic. This has never been 
the case here. However, it is clear that without the multiple-data 
per night policy adopted in the last periods of observations, it 
would be absolutely impossible to separate the rotations fourth 
harmonic from the alias of the period of CoRoT-7b. 
In what concerns the higher harmonics, we mention that the 
proximity of the period of the 6 th harmonic (3.94 days) to the 
period of CoRoT-7c (3.698 days) is a problem of difficult solu- 
tion. They are far enough to be separated one from another, but 
the beat period of the two components (about 60 days). is much 
larger than the timespan of the used subset (27 days). When the 
algorithm used in this paper is extended to include these har- 
monics, its convergence becomes excessively slow. The Fourier 
analysis of the residuals obtained with the 4-harmonic high-pass 
filter and an extended time interval (BJD 2,454,847 - 884) indi- 
cated that the amplitude of the radial velocities due to CoRoT-7c 
may be affected in up to 0.5 m.s^ 1 . This possible offset will be 
taken into account in the final results. 



4. The mass of the planets 

The technique described in the previous sections is certainly the 
best one we can devise to eliminate from a given series of un- 
evenly spaced observations the contributions of irregular long 
period terms and to get a remaining part which may be used to 



4 



S. Ferraz-Mello et al.: Planetary mass determination of super-Earths orbiting active stars 



determine the parameters of a short-period signal (the planets). 
It is now applied to the observations. However, for the reasons 
discussed above, we give up using all observations, but concen- 
trate on the set of observations done in the period BJD 2,454,847 
- 873. In addition to the aliasing problem, the consideration of 
the whole set - formed by 5 different subsets spanning 106 days 
plus 3 isolated observations one year before - is made difficult 
by the irregular variation of the stellar activity from the epoch of 
one set to the next. 

Another important setting in this determination is that we will 
concentrate on the masses and fix the periods in the values pre- 
viously determined (Leger et al. (2009) and QBM). One of the 
reasons is that having restricted the interval under study to a set 
of only 27 days, there is no possibility of improving a period 
determination resulting from observations taken from a time in- 
terval 4 times larger. The analysis of the covariance matrix with 
the 27- day data shows very high correlation (0.97) between pe- 
riods and phases for both planets what means that this short set 
cannot be used to simultaneously determine periods and phases. 
The results of this algorithm depend on the model used in the 
high-pass filter. Two main parameters were investigated. One of 
them is the model used in the interpolation to determine the ac- 
tivity at a given date: here, we considered both the 3- and the 4- 
rotation harmonics models. In a lesser extent, the 6-harmonics 
models has also been considered but the beat of the periods of 
CoRoT-7c and that of the 6 ,h rotation harmonic impairs the pro- 
cedure convergence. The other parameter is the Coherence Time, 
which sets the size of the window. We investigated several of 
them starting with the 20-day interval as used by QMB, but con- 
sidered also some other values in the range 8-22 days. The re- 
sults are shown in fig. [5] As far as the Coherence Time is kept in 
a limited interval, the results do not show large variations. Also, 
since the codes themselves depend on some operational param- 
eters that might affect the results, the procedures of filtering and 
mass determination were done with two very different codes: 
one, a lengthy steepest descent (diamonds) and the other a two- 
part code using a genetic algorithm completed with a downhill 
simplex (crosses). Labels indicate the corresponding Coherence 
Times. For the set shown by crosses, only the highest label (22) 
was shown to avoid excessive overlap of symbols and labels; the 
solutions with Coherence Times between 8 and 12 days cluster 
around Mb = 8raE art h. The others lie between these two limits. 
The green rectangle indicate a joint interval of confidence. The 
individual statistical errors of the self-consistent determinations 
were estimated as 0.5 m.s -1 that is, 0.7 and 1.1 Earth masses, for 
7b and 7c respectively. 

The sets of masses obtained for both planets show significant 
dependence on the used high-pass filter. The mass of CoRoT-7b 
shows some variations, following the 4 th rotation harmonic is in- 
cluded or not in the model. When the 4 f/! harmonic is included, 
we obtain for CoRoT-7b a mass above 8 Earth masses, while the 
results with only 3 harmonics is smaller than 8 Earth masses, On 
its turn, the models with 3 and 4 harmonics indicate, for CoRoT- 
7c, a mass around 14 Earth masses, but the results becomes 
less well determined when higher harmonics are included, in 
a way leading to believe that the actual mass of CoRoT-7c is 
smaller than the obtained value. The results of some runs using 
more observations (the 72 observations done in the interval BJD 
2454845 - 873), also result for CoRoT-7c a mass smaller than 13 
Earth masses. The results depend also on the period adopted for 
the planets. Because of its very short period, CoRoT-7b shows 
a larger sensitivity. Fortunately, the period of CoRoT-7b is very 
well known from the CoRoT photometry and this source of error 
can be discarded. 
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Fig. 5. Planet masses resulting from several runs of the iterative 
high-pass filtering with 3-, 4-harmonic filters (sets 3H, 4H). The 
other labels indicate the Coherence Time in days (see details in 
the text). The green square represents the interval of confidence 
of the solution given in Table 2. The large blue cross indicates 
the adopted solution in QBM and the red squares Ql and Q2 
the solutions corresponding to the published values of K\, K2, in 
QBM. The solutions A (alternative; only C0R0T 7b), B(Boisse), 
H(Hatzes), 5P (Fourier with 3 harmonics), FF(free Fourier) will 
be discussed in sections 8 and 9. 



The experiments done have also shown that the results are sensi- 
tive to the adopted weighting rules. All results in this paper were 
obtained using, at the beginni ng, the standard errors published 
by the observers (Leger et al. (l2009l) ). which were propagated 
following the classical rules of the least-squares formulas for 
unequally weighted observations (see Linnik, (|l963))Q. In the 
sequence, variances were obtained for each subset used in the 
filtering and were used to weight the averages on each date giv- 
ing the estimate of the activity. At last, let it be reminded that the 
mass determination from the filtered radial velocities cannot use 
the same weights as the given observations. The filtered RVs are 
differences between the observed RV and the estimated activity 
and are, therefore, affected by the errors in both these quantities. 
The classical formulas after which the variance (a 2 ) of the dif- 
ference is the sum of the variances of the two quantities entering 
in the subtraction is used and the new weights are defined as the 
inverse of the resulting variances. 



5. The star activity 

Figure [6] shows the activity determined using the 3-, 4-harmonic 
high-pass filters (bottom) and the difference between them (top).. 



4 The used classical weighting rules are able to take into account the 
fact that the observations done at the dates 2454860.75 and 2454864.63 
have a quality worse than the others; the weight associated to them was 
some 10 times less than the weight given to the more precise observa- 
tions in the set. 
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Fig. 6. Bottom: Star activity resulting from the high-pass filter. 
The labels [4H] and [3H] indicate the results obtained respec- 
tively with 4 and 3 harmonics. The color coding of the com- 
ponents of these curves, when visible, is the same as in Fig. [7j 
The dots indicate the measured radial velocities with their error 
bars. The error bars of the activity estimated with the 4 harmon- 
ics high-pass filter are also shown. Top: Difference between the 
results labeled 4H and 3H 
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Fig. 7. Top: Fourier Transforms (DCDFT) of the observed 
(black) and filtered data with the 3-and 4-harmonic high-pass 
filters in the period BJD 2,454,847 873. Bottom: Extension of 
the transforms to beyond the Nyquist frequency to show their 
behavior below 1 day. B, C are the orbital frequencies of CoRoT 
7b and 7c; B', C are aliases of B, C. The rotation period and its 
first harmonics are indicated on the top axis. 



The activities determined with both filters and various values of 
the Coherence Time are shown. For a given model, they do not 
show visible differences, notwithstanding the fact that these dif- 
ferences exist and affect the mass determination. The analysis of 
this figure may be summarized in a few words: (1) The filtering 
is very robust with respect to the chosen Coherence Time and 
model; (2) The results are more smooth when less harmonics 
are included in the filter. 

The activity estimated here may be compared to the one esti- 
mated by Pont et al. (1201 II) from the analysis of the bisector 
span measurements. It is remarkable that the features of the ac- 
tivity given if Fig. 2 of Pont et al. are very similar to those shown 
in fig. [6] and corresponding to the 4-harmonic filtering; however 
one may note that the total span of the RV due to the star activity 
is, there, about half of that shown in fig. [6] 

6. Quality of the new high-pass filter 

As done for the QBM determination, we may compute the 
Fourier transforms of the residuals V b s -V act i V and compare them 
to that of the observations. Some Fourier transforms (DCDFT) 
are shown in fig. [7] (top) for the two models and several val- 
ues of the Coherence Time. In order to avoid the complication 
arising from many almost overlapping curves, we present only 
those transforms for which the height of the peak near the fre- 
quency labeled B' (alias of the orbital frequency of CoRoT-7b) 
nearly matches the peak of the transform of the observed data. 
It is worth stressing that the frequency B' is not separated from 
the side lobes of the rotation period (which are broad because of 
the short time span of the used observation set). Because of this 
superposition, it is difficult to assess the quality of the filtering 



by inspection of the filtered spectra at this frequency. We may re- 
member that the composition of two frequencies in a spectrum is 
not just an addition since each of them carries one phase and the 
effect of the superposition cannot be assessed only by comparing 
their moduli. In this case, the superposition is reinforced by the 
short time span of the considered data. With a longer timespan, 
peaks would be sharper (as in fig. [2]) and could be separated, but 
using all available data in the analysis would mean to work with 
a discontinuous set of observations, which introduces additional 
(and in some extent unsolvable) difficulties in the estimation of 
the activity. One noteworthy effect of the superposition is the ap- 
parent enhancement of the peaks at B and B'. 
The frequency of CoRoT-7c, however, is less affected, at least as 
far as higher-order harmonics are not included. 
For the sake of giving an additional information on the aliasing 
effects, we present in fig. [7] (bottom) an extension of the trans- 
forms to an interval of frequencies including the actual orbital 
frequency of CoRoT-7b (labeled B) and one alias of the orbital 
frequency of CoRoT-7c (labeled C). As expected, the trans- 
forms in the given intervals are almost identical in both plots, but 
not perfectly identical because of the large proportion of nights 
with multiple observations in the selected time interval. 
The comparison of the Fourier transforms shows that the best 
quality filtering were obtained with Coherence Times 12-14 days 
when using the 4-harmonic filter and 18-20 days when using the 
3-harmonic filter (the curves for 18 and 20 days are almost iden- 
tical). 
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Fig. 8. Filtered radial velocities obtained with the 4-harmonic 
high-pass filter and Coherence Time t = 14d (top) compared to 
those given in QBM (bottom). 



7. The filtered radial velocities 

Figure [8] (Top) shows the radial velocities obtained subtracting 
from the measured radial velocities the activity resulting from 
the use of the 4-harmonic filter and Coherence Time r = 14d. 
They may be compared to the ones used in the discovery paper 
(Fig. [8] Bottom). The main differences appear in the beginning 
of the interval, where the filtered radial velocities appear much 
larger in our results than in QBM and in the middle of the inter- 
val where the contrary occurs. It can be easily seen that around 
the date 42, in our results, a destructive beat occurs between the 
RV sine curves corresponding to the two planets, at the actual 
observing times. 

We should stress the fact that the epoch BJD 2,454,820 used in 
the plots and calculations of this paper is not the same epoch 
used in other papers. It is an arbitrary epoch close to the ac- 
tual dates and allows the phases (i.e. the longitudes at the fixed 
epoch) to depend less strongly on the periods, thus making the 
numerical procedures more robust. In adition, we note that in 
all steps one additive constant representing the radial velocity of 
the planetary system is determined together with the other un- 
knowns. 



8. A model-independent approach 

The use of a high-pass filter allowed us to separate parts of 
the RV measurements due to high and low frequency compo- 
nents. The results are, however, model-dependent. The high- 
pass filtering has given statistically coherent estimates for the 
mass of CoRoT-7c, but that result can carry some systematic 
effects due to some higher-order rotation harmonics. In what 
concerns CoRoT-7b, the filtering led to two different solutions. 
From the purely statistical point of view, the 4-harmonic filter- 
ing should be better than the 3-harmonic filtering. However be- 
cause of the aliasing involving the rotation A' h harmonic and the 
orbital period of CoRoT-7b, it is convenient to confirm the re- 



sults with some alternative technique independent of assump- 
tions involving the rotation period and its harmonics. An al- 
ternative approach was suggested by Hatzes et al. (2010) (and 
paper in preparation) in which only data from nights in which 
multiple observations were done are used, and assuming that in 
the 4 hours time span of one observation night, the activity of 
the star does not change and may be fixed as the same for the 
three observations . We note that there are 10 nights in the in- 
terval 2,454,847 - 873, in which three consecutive observations 
were done. Then, the problem is to fit those observations with 2 
planetary sine waves (the eccentricities may be assumed equal 
to zero as discussed in section 10) and a set of 10 independent 
constants added to the data in the corresponding dates. The main 
problem with this approach is the small number of degrees of 
freedom in the best-fit problem. We have 30 data and 14 un- 
knowns. The number of degrees of freedom is 16. In analogy 
with a chi-square distribution, we may guess that every solution 
leading to a w.r.m.s <0.6 m.s higher than the minimum shall 
be considered as belonging to the standard interval of confidence 
(see Press et al. d!986l) ). In order to improve the determination, 
we have chosen to add to these observations those from 16 other 
nights in which 2 observations were done. We have then 62 data 
and 30 unknowns. The number of degrees of freedom becomes 
32. This is twice of what we would have if only using the 10 
nights with 3 observations each and means a less broad confi- 
dence interval, including only the solutions leading to a w.r.m.s. 
<0. 16 m.s -1 larger than the minimum. 

To assess the set of all these 'good-fit' solutions, we may use 
the same biased Monte-Carlo technique used to obtain good fits 
for the planets HD 82943 b,c (Ferraz-Mello et al. d2005l) ). As be- 
fore, the frequencies were fixed at the values given in QBM (they 
are well determined) and only 4 unknown planetary parameters 
were considered: The two masses and the two phases. The (bi- 
ased) random search produced thousands of solutions. Those 
with w.r.m.s. < 1.67 m.s 1 are shown in fig. [9] In that figure 
we superpose the results for the half-amplitudes K of the two 
planets. The mass of CoRoT-7b is constrained to the interval 
niB - 7.2 + l.4mEanh (i.e.4.l < K < 6.2 m.s~') and its longi- 
tude at the epoch BJD 2,454,820.0 is 200 + 7 deg. The mass of 
CoRoT-7c, however, is not constrained by this set of observa- 
tions. Fits with w.r.m.s. < l .67 m.s~' exist for masses in a broad 
interval: mc sin ; < 22mEarth (Kc < 9.8 m.s~'). 
The poor constraining of the solution for CoRoT-7c deserves 
some discussion. The alternative technique proposed by Hatzes 
et al. d20l0l) is good to separate low and high-frequency contri- 
butions. The period of CoRoT-7c is not short enough to be con- 
sidered as low. In the more favorable situations, the variations 
of the RV due to CoRoT-7c may reach 2 m.s -1 in 5 hours. It is 
generally less, and not distinguishable from contributions com- 
ing from low-amplitude long-period terms. We cannot avoid the 
random process used from mimicking, in the activity, contribu- 
tions with the same period as CoRoT-7c and from combining 
them with the planet RV. 

The results are not accurate but they are important because of 
their independence with respect to assumptions on the activity 
behavior. 

The best constrained fit, defined as the leftmost point in fig. [9] 
corresponds to the mass mj = 7.5 Earth masses. This fit and the 
interval of confidence corresponding to the maximum span of 



5 The results from the used RV modeling indicate that variations are 
smaller than the observational errors. However, in one extreme (and 
infrequent) case, the estimated activity varied by 0.8 m.s~' between the 
first and the last observation done in the night. 
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Fig. 9. Good fits obtained with a biased Monte-Carlo under the 
condition wrms < 1.67 m.s -1 . The blue dots correspond to 
CoRoT-7b and the red ones to CoRoT-7c. 



mthe solutions oin the considered interval is shown by point A 
in fig. H 

It is important to stress that each set of masses and phases of the 
two planets obtained in this experiment corresponds to one solu- 
tion whose residuals fit the observations better than 1.67 m.s -1 
(w.r.m.s.) and cannot be discarded. The density of the points in 
fig.|9]has no meaning (it rather reflects the strategy used to con- 
struct the sets). We note that only the error bar for ms is given. 
It was not possible to get a good determination for m c sin i (the 
error bar would extend more or less over the whole width of the 
figure). 

9. Fourier approaches 

As a first guess, one could assume that the observations can 
be represented by the sum of several periodic functions and 
use conventional non-windowed Fourier analyses to determine 
their parameters (amplitudes, periods, phases). The solutions 
ind icated by B and H in fig. [5] are those obtained by Boisse et 
al. (1201 ll) and Hatzes et al. y010) from the Fourier analysis of 
the observations. Boisse et al's solution is the best fit of 5 sine 
curves to the radial velocities, 3 of them having the periods of 
the rotation and its first harmonics. For that solution, besides 
the 52 points used in this paper, the remaining 9 points of the 
latest observational subset were also considered. Hatzes et al's 
solution was obtained using all observations and a pre-whitening 
plus filtering procedure based on periodograms; 9 different sine 

curves contributed to the construction of the solution. 

If we proceed in the same way as Boisse et al. d2011l) but 
considering only the set of 52 observations done used in 
previous sections, the result is that shown by the symbol 
5P in fig. [5] In fact, we have done a great deal of Fourier 
analyses of the measurements, which remained unpublished 
only because they assume a periodic behavior in the activity 
which is very improbable. Every function in a finite interval 
may be represented by a Fourier series, but in the present case, 
the main period needs to be the rotation period. Such periodic 
structure is not seen in the residuals published by QBM (shown 
in fig. Q] above), in our own activity curve (shown in fig. [6] 
above), in the activity as determined by Pont et al.( (1201 ll) : 
their fig. 1) or in the photometric series produced by CoRoT 
(Leger et al. (2009)). In all these plots there is some kind 



Table 1. Fourier decomposition of the data in the time interval 
2,454,847 - 873. The resulting w.r.m.s of the residuals is 1.76 
m.s -1 



Period (d) 


K(m.s-') 


Mass (Earth=l) 






x sin i 


21.16 + 0.6 


15.0™ 




3.70 ± 0.04 


6.0 ± 0.9 


13.4+ 2.0 


0.850 ± 0.002 


5.8 ± 1.2 


8.1 ± 1.5 


1 1 7 +0 3 

1J --' -0.7 


4.0 ± 1.0 




5.09 ±0.15 


2 4 +lA 





of repetition associated with the rotation period, but in some 
sections the curve appears dominated by the rotation period, in 
others by the rotation harmonics, and so on. None of the partial 
curves appearing in fig. Q] shows the 23.64 period indicated 
by photometry. An irregular behavior is also seen in the light 
curve, whose periodogram clearly shows the rotation period and 
its harmonics, but which is such that in some sections of it no 
variations are seen, while in others the rotation period is well 
marked. 

We report here only one of the experiments which consisted 
of several steps: (a) The periods were initialized in relatively 
broad intervals bracketing the rotation period and its harmonics, 
allowing for possible differences due to either the physics of the 
activity or the beat between the actual periodic signal and the 
sampling set (actual observation dates), (b) The approximated 
more probable elements (including the periods) were determined 
via a chain of best-fits to the data; (c) In each step, the elements 
were determined by the simultaneous best-fit of N trigonometric 
components. The process was stopped when the addition of one 
new term was no longer able to improve significantly the results 
(F-test) and the scrambled data produced spectra with peaks of 
the same size as the ones obtained with the real data. In such 
case, the inclusion of more unknowns in the process may lead to 
high correlation between the unknowns and to results that may 
become undistinguishable from artifacts. 

The eventual results of the free Fourier approach are shown 
in fig. [5] (point FF) and detailed in Table Q] The error bars 
were estimated with a biased Monte Carlo sampling of the 
neighborhood of the solution, which showed very asymmetric 
distributions in some cases. 

It is important to stress the fact that the weakest of the five 
terms is already very uncertain as both the spectrum of the 
residuals in the previous step and the a posteriori F-test have 
shown that no actual improvement was obtained by the addition 
of this term to the solution. 

One important by-product of this analysis concerns the con- 
jecture of the existence of one third planet at P — 9.2d, which 
cannot be studied with the other techniques described in this 
paper, since they do not allow such a slow periodic variation 
to be distinguished from the star activity. Peaks corresponding 
to periods around 9 days are recurrent in all analyses done, 
since the beginning of this investigation. They can be seen 
for instance in the spectra of the raw data in fig. [2] and fig. [7] 
However, always, when a monochromatic filter is used and a 
sinus curve with the amplitude and period of the rotation is 
subtracted, that peak disappears. We conjecture that it results 
from a complex beat between the rotation period, the sampling 
dates and the planets themselves. 
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Fig. 10. Orbit evolution of two planets with masses 8 and 15 
Earth masses in a CoRoT-7-like system with mutually inclined 
orbits. The damping of the eccentricities is due to the tides on 
the inner planet only. 

10. The eccentricities 

In the mass/orbit determination presented in this paper, the ec- 
centricities were taken equal to zero. Indeed, an analysis of the 
dynamical problem shows that the tidal dissipation in the planet 
CoRoT-7b and its gravitational interaction with CoRoT-7c damp 
both eccentricities. Simulations were done where t he forces due 
to the tide raised in the planet (cf Mignard dl979l) ^ were added 
to the gravitational ones. In all of them, the osculating eccentric- 
ities of the two planets stabilize in values of the order of resp. 
1CT 5 and 1CT 4 in a few tens of Myr, whichever initial eccen- 
tricities and inclinations are considered. In the solution shown in 
fig-HQ] the planets are initially on circular orbits with semi-major 
axes 0.0175 and 0.0456 AU (like CoRoT-7c) in two planes with 
a mutual inclination of 40 degrees. The chosen initial eccentric- 
ities do not affect the solution because, immediately after the 
beginning of the simulation, the gravitational interaction of the 
2 planets forces the eccentricities to be larger than 0.1 and 0.01 
respectively (In the coplanar simulations they jump to the equi- 
librium values fixed by the mutual perturb ations" cf. Mardling 
d2007l) : Rodriguez d2010l) ; Rodriguez et al. d2011l) ); Thereafter, 
the tidal dissipation in the inner planet starts dissipating the en- 
ergy of the system, making the orbit of CoRoT-7b slowly spiral 
down towards the planet and become circular while the semi- 
major axis of CoRoT-7c remains almost unchanged. In addition, 
the exchange of angular momentum between the two planets also 
damps the eccentricity of CoRoT-7c and drives the planets to an 
equilibrium configuration (see Mardling (1201 ll) ~). It is worth re- 
calling that after the circularization of the orbit of CoRoT-7b, the 
tidal f riction in this planet almost ceases (see Ferraz-Mello et al. 
d2008ln and no longer continues to significantly affect the evolu- 
tion of the system. Let be added that in some coplanar runs start- 
ing with eccentric orbits, the final eccentricities are yet smaller 
than the ones shown in fig. ITUI 

11. Discussion 

Three independent approaches were used to assess the masses 
of the exoplanets CoRoT-7b and CoRoT-7c notwithstanding the 
difficulties created by the activity of the star which contaminates 
the measured velocities not as a jitter, but as a dominant signal 
some 2-3 times more intense than the planetary contributions. 
Before comparing the results, it is worth stressing the fact that 




30 40 50 
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Fig. 11. Bottom: Filtered radial velocities with the 4-harmonic 
high-pass filter and 14-day Coherence Time (dots) and RV curve 
corresponding to the planets solution obtained with the same pa- 
rameters. Top: Differences between the filtered RV and the RV 
curve. 



the residuals of the 3 analyzed models have weighted r.m.s. less 
than 1.9 m.s~', what may be considered good given the irregular 
activity of the star. However, we have to prevent against using 
the minimum w.r.m.s. of the models to mutually compare 
them, because, by definition, these quantities are minima of 
the residuals with respect to a given model and therefore, may 
be considered as model-dependent. They shall be used only to 
make comparisons within a given model. 

The introduction of a self-consistent algorithm improved the 
filtering used in QBM and resulted in masses considerably 
larger than the ones previously obtained. One fact influencing 
the results is the adoption of a 4-harmonic filter in addition to the 
3- harmonic filter. The analysis of CoRoT's photometric series 
did not reveal the need of using the fourth harmonic. However, 
the analysis of the spectrographic parameters associated with 
the activity of the star (mainly log 7?^) showed the importance 
of high-orders harmonics even beyond the fourth. In addition, 
the fact that the 4 th harmonic is very close to an alias of the 
period of CoRoT-7b, makes its introduction necessary. 
The coincidence between the period of the fourth rotation 
harmonic (5.91 d) and an alias of the orbital period of CoRoT-7b 
(5.925 d) is a problem of major concern in this study. We have 
to stress first that in series of astronomical observations done 
at the same hour angle it is impossible to solve any aliasing 
problem. In such case, a linear relationship appears between 
the two components. We may force at will the value of one 
component and compensate with the value of the other. There 
is no mathematical tool able to solve this problem. In this case, 
given the slight difference between the two frequencies (5.91 
and 5.925) and the complexity of the signal, the only way to 
get one solution is to break the uniformity of the time intervals 
in the sample by means of observations done in very different 
hour angles. Unfortunately, with observation from only one 
observatory this is not easy. Differences of a few hours with 
respect to the mean are the maximum that could have been 
obtained in this case. The set of observations selected for this 
study has 52 observations done in 27 nights, including 10 nights 
with 3 observations in a 4 hours interval and 5 nights with 2 
observations. This is the only hope that we have to solve the 
beat of the two frequencies before new observations are done. 
They were indeed separated using a 4-harmonic high-pass filter 
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and the convergence of the self-consistent iteration routes to 
the same result (see fig. [3) is an indication that, in the adopted 
algorithm, they may be considered as independent. 
One additional comment with respect to the beat comes from the 
analysis of fig. |6](top). That figure shows the difference between 
the estimated activity in the two cases: with 3 and 4 harmonics, 
respectively. It is an evaluation of the contribution of the fourth 
harmonic to the estimated activity of the star. The surprise, in 
this case, is that the apparent period of this contribution is not 
5.91 d, but only 5.25 d. Assuming that the rotation period is 
the same observed by CoRoT, this difference would mean that 
we are not dealing with one frozen periodic signal, but with a 
signal that is just nearly periodic: as the evolution of a periodic 
process whose period and phase are continuously changing. 
The other possibility is that the period is not the same as the 
published one. We note that 5.25 is one fourth of the main 
period found in the free Fourier analysis of the radial velocities, 
21.15 days, and that the frequency of the highest peak in the 
Fourier spectrum of the data, in fig.[7]is near 0.5 d F -l instead of 
0.042 d -1 , inverse of 23.64 d0 Whatever the reason responsible 
for this difference, it certainly contributed to the fact that the 
4-harmonic filter succeeded to get one solution with the two 
components separated, while classical least-squares solution 
and a covariance analysis using the filtered RV obtained with the 
published rotation was unable to separate the orbital period of 
CoRoT-7b from the alias of 4 f/! harmonic. As a check, we have 
done some runs of our codes using this lower period. The results 
are masses near the values obtained using the actual values of 
the rotation period. So they confirm those results and show that 
the filtering algorithm is robust with respect to variations in the 
rotation period used in the filter, at least as far as it is close to 
the adopted ones. 

The quality of the results may be assessed from the residuals 
(O-C) and from the fitting of the observations to the model 
shown in fig. QT| The obtained (O-C) correspond to a w.r.m.s. 
1.9 m.s -1 . Two comments may be added here. (1) The error 
bars in fig. QT| are larger than those of the RV measurements. 
They result from the addition of the variances of the RV 
measurements to that of the activity given by the used model. 
The large errors in some dates in the middle of the interval 
comes from the fact that, in those dates, only one observation 
was done per night, leading to larger statistical errors in the 
estimation of the activity in these dates; (2) These results absorb 
naturally RV jumps as those that led Pont et al (2010) to state 
that HARPS systematic errors may be huge in some cases. 
One of the 20 m.s 1 jumps reported by them, which occurred 
between the dates BJD 2454868/69 (BJD - 2454820 = 48/49), 
is in the studied interval. However, in our results, it is mainly 
the consequence of an important increase in the star RV due to 
the reaction to the combined motion of the 2 planets at those 
dates. When the increase due to the activity (about 3 m.s -1 ) is 
added to it, one may verify that the part of the RV increase non 
explained is still large, but just 1/3 of the increase seen in the 
RV measurements. 

An alternative approach suggested by Hatzes et al. (1201 Oh . 
which has the advantage of making no explicit hypotheses on 
the activity (except its slow variation) has also been improved 



6 The variation from 23.64 to 21,1 days in the activity period is con- 
sistent with differential rotation differences expected in a solar-type star 
(see Thomas and Weiss, i2008l) ). However, the time span of the observa- 
tions here considered is not large enough to allow us to give full credit 
to this result and the question should be reconsidered when new data 
become available. 
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Fig. 12. Longitudes at the epoch B JD=2,454,820 resulting from 
several runs of the iterative filtering. The labels indicate the num- 
ber of harmonics in the high-pass filter and the Coherence Time 
in days. For the results obtained using the genetic algorithm + 
simplex code (crosses), only some results are labeled to avoid 
overlaps. 

Table 2. Masses and elements (Epoch BJD 2,454,820.0) 





CoRot-7b 


CoRot-7c 


Kfm.s- 1 ) 


5.7 ±0.8 


6.0 ±0.6 


Mass (Earth Mass) (t) 


8 ± 1.2 


13.6 ± 1.4 


Longitude at Epoch (deg) 


207 ±4 


225 ±3 


Period (d) (fixed) 


0.853585 


3.698 


Period (d) (estimated) 


0.85354 (see text) 




TO (BJD) (estimated) 


2,454,847.893 ± 0.009 




Eccentricity (fixed) 









t The mass of CoRoT 7c is minimal (x sin i) 



and combined with a biased Monte Carlo approach. The 
resulting mass of CoRoT-7b lies close to results obtained using 
the 3-harmonic filter. The resulting mass of CoRoT-7c is ill 
determined because the period of CoRoT-7c is not small enough 
to allow it to be fully separated from the star activity. 



12. Conclusions 

The final results were obtained by combining the results of the 
3- and 4-harmonic filtering and of the multiple-observations-per- 
night alternative approach. They are given in Table [2] The re- 
sults from filtering were considered when the two codes used 
were convergent and lead to nearby values. Otherwise they were 
not considered in the final estimation of the results. In particu- 
lar we mention that the largest and least Coherence Times (22 
and 8 respectively) were not considered. The other criterion was 
the quality of the filtering as indicated by the comparison of the 
Fourier Transforms before and after the filtering (fig.|7) 
We remind that, because of the small timespan of the obser- 
vations used, we renounced to determine the periods and used 
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those given in the discovery papers (Leger et al. (2009]) and 
QBM). 

With the period determined from the transits and the date of the 
first transit of CoRoT-7b observed by CoRoT, we obtain 196.4 
degrees for the expected longitude at the epoch BJD 2454820.0. 
The small offset with respect to the result shown in tabl43]corre- 
sponds to a slight correction in the period, which would then be 
0.85354 + 0.00002 days. 

We also stress the fact that the errors given in table [2] intend 
to define the interval in which the two masses are expected to 
be. They are not the result of some popular statistical formu- 
las. Those formulas generally are valid in conditions much more 
restrictive than the actual ones and often the results are much 
smaller than the actual errors. An additional difficulty in the re- 
sults for CoRoT-7c comes from the fact that it is difficult to dis- 
entangle the mass of this planet from the rotation higher har- 
monics (mainly the 6 th ). In the case of the beats between the 
period of CoRoT-7b and the 4 th harmonic, it was possible to dis- 
entangle them by making multiple observations in many of the 
nights of the latest observation periods. In the case of CoRoT-7c, 
the beats are not related to aliasing, and the problems they raise 
cannot be solved by the strategy of making multiple observa- 
tions per night nor by mathematical tools allowing low and high- 
frequency components to be separated. To solve this kind of en- 
tanglement we should have observations spanning over a larger 
time interval. We have done some calculations using all observa- 
tions available in the interval BJD 2554825 - 884, increasing the 
time span from 27 to 60 nights. This means to include discon- 
tinuities in the observations set, which make the analysis more 
difficult and, in some sense, of less confidence. In these runs the 
mass obtained for CoRoT-7c is below 13 Earth masses. Adding 
this to the evidences from the Fourier analysis of the residuals 
used in the runs labeled as 4H, that the 6 th rotation harmonic 
may contribute with 0.5 m.s -1 to K2, we may guess that values 
between 13 and 14 Earth masses are the more probable ones. 
One may wish that the improvement of techniques as the one 
recently proposed by Pont et al. (1201 ll) applied to simultaneous 
photometric and spectrographic observation done over long time 
spans be able to give a final answer to this question in future. 
However, with the existing observations, we are confident that 
the above estimates are the best ones we can obtain and that all 
consistent estimations fall in the ranges given above or, at least, 
in its immediate neighborhood. 

The longitudes at BJD=2,454,820 were obtained in the same 
way as the masses, comparing the solutions from various mod- 
els (see fig.fT2l. The individual statistical erorrs are 5 and 12 de- 
grees for planets C0R0T 7b and 7c respectively. The other two 
approaches did not contribute to the given results because of too 
large errors. 

One last point to note is that the result found here for C0R0T 7b 
means that the planet has a bulk density 11 +3.5 g/cm 3 . This is 
much more than the 6.6 g/cm 3 resulting from the mass given in 
QBM and mean not only that C0R0T 7b is rocky, but also that 
the contribution due to the iron core must be higher (between 50 
and 6 5 percent in mass; see Fortney et al. (2007); Seager et al. 
( 2007 )) and that the density at its center may be close to 25 g / cm 3 
(see Seager et al. (120071) 1. These values are higher than those 
corresponding to other known rocky planets, but all consistent 
determinations of the mass of CoRoT-7b lead to bulk densities 
of at least 9g/cm 3 showing that the high-density of CoRoT-7b 
is a constraint to be taken into account in the modeling of the 
planet. 

The comparison of the various methods used in this paper al- 
lows us to say that the high-pass filter used by QBM (Queloz et 



al. d2009)), embedded in a self-consistent algorithm, is the best 
one we can devise to disentangle long- and short-period terms 
in a given series of unevenly spaced observations. Its superiority 
over Fourier analyses with fixed frequencies comes from the fact 
that the used running window allows the method to treat a signal 
which is not periodic or has a period different of the pe riod u sed 
in the filter. The procedure proposed by Hatzes et al. (1201 Oh of 
using only nights with multiple observations and including an 
additional Vq for each date is an important complement. It suf- 
fers, however, from a shortcoming due to the great number of un- 
knowns that it involves. As a consequence the number of degrees 
of freedom is small and the resulting confidence intervals are too 
large. In the current case, it just allowed us to determine the mass 
of C0R0T 7b. This was great because this was the ultimate goal 
of our work, but the impossibility of getting reasonable confi- 
dence intervals for the mass of C0R0T 7c shall be mentioned. 
The planning of new observations should take this into account 
and have as many multiple night observations as possible. 
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